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Abstract 

We construct a class of quantum stochastic models of reservoir driven many-particle 
systems that are the natural counterparts of certain extensively studied classical ones, 
which have been shown to exhibit good hydrodynamical behaviour. Our treatment of these 
models achieves two main aims. The first is to show that they enjoy the hydrodynamical 
properties of their classical counterparts. The second is to show that they satisfy the key 
assumptions of the general quantum macrostatistical scheme, presented in earlier works by 
the author, which served to expose certain generic large scale features of nonequilibrium 
steady states, e.g. the long range hydrodynamical correlations that they carry. In this way 
we establish the viability of that scheme. 
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I. Introduction 



The purpose of this note is to bring together two developments in the theory of the 
relationship between the hydrodynamical and the microscopic pictures of reservoir driven 
macroscopic systems.. The first of these developments, which we shall refer to as (I), is a 
body of work concerned with the derivation of hydrodynamics from the microscopic dynam- 
ics of a class of classical stochastic models [1-5]. Most interestingly, the hydrodynamical 
fluctuations of these models about their nonequilibrium steady states have been shown to 
carry long range spatial correlations [3-5] and to conform to a generalised version [5] of the 
Onsager-Machlup process [6]. The second development, which we shall refer to as (II), is 
a general, model-independent, quantum macrostatistical treatment [7,8] of hydrodynami- 
cal fluctuations about nonequilibrium steady states, that is based on certain hypotheses of 
chaoticity and local equilibrium, together with a generalised version of Onsager's regression 
hypothesis [9]. On this very general basis, we have shown that, as in the special classical 
models of [3-5], the hydrodynamical fluctuations about nonequilibrium steady states carry 
long range spatial correlations and execute a generalised Onsager-Machlup process. 

Our objectives here are to extend the constructive developments of (I) to the quan- 
tum regime and to show that the resultant models satisfy the assumptions of the general 
quantum macrostatistical theory of (II) . The latter objective is thus designed to show that 
the 'axiomatic' scheme of (II) is viable. 

Our approach to these objectives is based on a construction whereby we extend the 
generic classical stochastic model, E c ;, of (I) to a quantum system, E, in such a way that 

(i) the abelian algebra of observables, £>, of E c ; is a subalgebra of the nonabelian one, A, 
ofE; 

(ii) B is stable under the dynamics of E; 

(iii) the nonequilibrium steady state of E c ; is just the restriction to B of that of E; and 

(iv) the hydrodynamical observables of E are precisely those of E c ;, 

Thus, the construction of the quantum model E in this way permits us to exploit some of 
the powerful results obtained for the hydrodynamical properties of its classical counterpart, 
E c j. In particular, it enables us to verify that E enjoys the hydrodynamical properties of 
E c j and, moreover, that it satisfies the basic assumptions of (II). This latter result therefore 
establishes that the scheme (II) is viable. 

We present the formulation of the models E c j and E in Section 2, and establish there 
the above properties (i)-(iv). Thus the hydrodynamical picture of E reduces to that of 
E c j. We formulate this picture explicitly in Section 3 and extend it to fluctuations about 
a nonequilibrium steady state in Section 4, where we specify the regression, chaoticity and 
local equilibrium hypotheses on which the theory of (II) was based. In Section 5 we prove 
the validity of these hypotheses for the present stochastic model and thereby establish the 
viability of the scheme of [II]. We conclude in Section 6 with a brief comment about an 
outstanding problem in the theory of quantum stochastic processes. 
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2. The Classical and Quantum Stochastic Models 



The generic model with which we shall be concerned, whether classical or quantum, is a 
system of N identical particles that live in a bounded region Ojv of the ^-dimensional lattice 
Z d and are coupled to reservoirs at its boundary. More specifically, O^v is assumed to be the 
subset of Z d contained within the dilation by a certain factor, L^v, of a fixed, iV-independent 
bounded open connected region O of the Euclidean space R d . Thus, 0_?v = Z d n(LjvO). We 
define its boundary, <90jv, to comprise the sites in O^r with at least one nearest neighbour 
that lies outside that region, and we define Int(Ojv), the interior of O N: to be Qn\9^n- 
This latter region thus consists of the sites in fl N whose nearest neighbours also lie in Q N . 

We assume that the volume of O is unity, that its boundary, <90, is smooth and that 
the mean particle number density, z/, of the system is iV-independent. Thus 



We assume that the dynamics of the model corresponds to a stochastic process whereby 
the particles jump between nearest neighbouring lattice sites according to probabilistic 
laws that will be prescribed below. 



For the classical model, E c ;, we denote by n x the number of particles at the site x. 
In the case where an exclusion principle is operative, n x is restricted to the values and 
1: otherwise it may take any non- negative integral value. Thus a particle configuration is 
a map n : x^>n x of O^ into a set K, which is either {0, 1} or N according to whether or 
not an exclusion principle is operative. We take the algebra, £>, of bounded observables 
of the system to comprise the bounded, complex valued functions on the configuration 
space T = K N , with supremum norm. Thus, equipping Y with the discrete topology, 
B = C(r), the C*-algebra of bounded continuous functions on T. For x, yefijv, we define 
n-^n x ' y to be the transformation of V corresponding to the transfer of a particle from 
x to y, provided that that transfer is kinematically admissible, i.e. that (n x — 1) and 
(n y + 1) lie in K: otherwise we define n x,y to be simply n. Likewise we define n x ^ to be 
the modifications of n corresponding to increments ±1 in n x , provided that (n x ±l)&K: 
otherwise we define n x,± = n. We assume that the dynamics of the system is given 
by a continuous one-parameter semigroup, c j(R+) = {0 c /(£)|teR+} of linear, positivity 
preserving transformations of B. We denote its generator by Q c \ and we shall presently 
specify its explicit form for two models, namely those conventionally termed [1-5] the 
simple exclusion model and the zero range model. 

The Simple Exclusion Model. . For this model, K = {0, 1} and Q c \ takes the following 

form 



where the prime over the first sum signifies that summation is confined to nearest neigh- 
bours, h is a smooth positive-valued function on 90 and r& is the number of nearest 



L N = (N/u) 1 ^. 



(2.1) 



2.1. The Classical Model. 
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neighbours of 6(g90at) on the lattice Z d that lie outside Ojv- Thus, the first sum repre- 
sents the jumps between nearest neighbouring sites of the particles in the interior of Qn, 
the second the escape of particles across its boundary and the third the supply of particles 
by external sources at the boundary. 

The Zero Range Model. For this model, K = N and, in the same notation as in Eq. 
(2.2), Q c i takes the following form. 

Qdf(n) = Y! PO 9M{f(n x ' y )-f(n)) + 

£ 6eaiJ /* M/w - /(*)) + E 6eafi „ WLiv) (/(<) - (2 - 3) 

where g is a positive valued, non-increasing function on K for which g(0) = and 
sup fc ((7(/c + 1) — g(k)) is finite. 

Note. It follows from Eqs. (2.2) and (2.3) that in the cases of the simple exclusion 
and the zero range models Qn x takes the forms (An)j for the and (A^on))^, respectively, 
for xelnt(Ojv)) where A is the discrete Laplacian defined by the formula 



< * -7/ f= S I AT 



ye^N 

and the prime over S again indicates that the sum is taken over sites y that are the nearest 
neighbours of x. Hence, for both models, the dynamics of the field n is diffusive. 



2.2. The Quantum Model 

We take the quantum model E to be a system of fermions or bosons according to 
whether or not the exclusion principle is operative. In either case we formulate the model 
in a standard way in terms of the Fock space H and the creation and destruction operators 
{a*, a x |a;GOi} that act therein according to the following defining conditions. 

(a) TC contains a vector $ that is annihilated by the action of each of the a x 's and is cyclic 
with respect to the polynomials in the a*'s; and 

(b) the operators a x and a* satisfy the canonical commutation or anticommutation rela- 
tions, namely 

[ax,a*] T = 5 x , y I; [a x ,a y ] T = V x, yeVL N: (2.4) 

according to whether the system consists of bosons or fermions. For either case, we define 
the number operator 

h x = a*a x V xgOat. (2-5) 

It follows immediately from Eqs. (2.4) and (2.5) that the h x s intercommute and thus con- 
stitute a classical field n := {ti x \xeQn}- We denote by ip(n) the simultaneous eigenvector 
of these operators h x with corresponding eigenvalues n x , i.e. 

n x ip(n) = n x ip(n) V xE^In- (2-6) 
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It follows from this formula and our specifications of H that the vectors ip(n) form a 
complete orthogonal basis for this space as n runs through the classical configuration 
space r = K N , with K = N or {0, 1} according to whether the particles are bosons or 
fermions. Further, by Eqs. (2.4)-(2.6), 

a x %l){n) = nl /2 ij(n x >-) and a%^{n) = (l±n x ) 1/2 i(j(n x > + ) (2.7) 

and hence 

a^aMn) = (n x (l±n v )) 1/2 i(>(n x >v), (2.8) 

where n x,y and n x,± are defined as in Section 2.1 and ± signifies the boson-fermion alter- 
natives. 

We denote by T the additive group of bounded continuous real-valued functions 6 : 
x^9 x on Oat and we define the unitary representation U of T by the formula 

C/(0) = exp(zV 9 x n x ), (2.9) 

We take the algebra, A, of bounded observables of E to be that of the bounded 
operators in H and we define 7 to be the representation of T implemented by U in Aut(*4.), 
i.e. 

-y(9)A = U{e)AU{e)~ 1 V AeA. (2.10) 

Thus, by Eqs. (2.4)-(2.6), (2,9) and (2.10), 7(6*) is the local gauge automorphism given by 
the formula 

l(9)a x = a x exp(-i9 x ) (2.11) 

We define B to be the locally gauge invariant subalgebra of A i.e., by Eq. (2.10), the 
set of elements of A that commute with all the L r (6') , s. It follows from this definition and 
Eq. (2.9) that B comprises the elements B of A for which 

U(6)Bip(n) = BU(9U(n) = exp(zV n x x )Bi;(n) V neT, OeF, 

i.e. for which Btfj(n) is a simultaneous eigenvector of the h x s with corresponding eigenval- 
ues n x . This signifies that Bip(n) = F(n)(/)(n), where F is some bounded complex-valued 
function on T, i.e. that B is the operator F(h), defined by the formula 

F(n)ip(n) = F(n)ip(n) V neT. (2.12) 

Thus we have established the following proposition. 

Proposition 2.1. (1) B comprises the functions of n; and 

(2) the mapping F(h)—>F(n) of B onto C(T), the algebra of observables of 'Ed, is a C* - 
isomorphism. Hence B may be identified with the latter algebra. 
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We assume that the dynamics of the model £ is given by a strongly continuous one- 
parameter semigroup 0(R+) = {(p(t)\t£Tt + } of completely positive identity preserving 
contractions of A. Its generator Q therefore takes the standard form for that of a quantum 
dynamical semigroup, namely [10, 11] 

QA = i [H, A] _ + ^{VfAVj - \\VfV Si A]+) V AeA, (2.13) 

where H is a self-adjoint element of A and Vj and "Y^jV^Vj a ^ so belong to this algebra. 
Thus, the dynamics of the model is determined by H and the V^'s. We shall now specify 
thses operators for quantum versions of the simple exclusion and zero range models. 

The Quantum Simple Exclusion Model. In view of the requirement of an exclusion 
principle, we take the particles of this model to be fermions. We construct its dynamical 
semigroup 0(R+) in such a way as to obtain a natural correspondence between its generator 
and that of C /(R+), as given by Eq. (2.2). Specifically we choose the operators H and 
the V : s of Eq. (2.13) in the following way. 

(a) We take H to be zero, since Q c \ contains no Hamiltonian part. 

(b) Since the indices involved in the structure of Q c \ comprise the nearest neighbouring 
pairs (x, y) of sites of Qn together with the boundary sites b of <9Qat, we assume that the 
index j of Eq. (2.13) also runs through just these sets. 

(c) For j = (x, y), we choose Vj to be a*a x , since the first summand of Eq. (2.13) then 
represents the transfer of a particle from x to y, with probability rate that corresponds to 
that of Eq. (2.2). 

(d) For each &G<9fi, we introduce two separate V's, namely h{b/ L^) 1 / 2 ^ and r\j 2 a\,, 
which lead to the creation and annihilation, respectively, of the particle at 6, with weights 
corresponding to those of Eq. (2.2). 

Thus, under these specifications, the formula (2.13) takes the following form. 

gA = Y] (a* x ayAa* y a x - ha*a y ala x , A}+) + 

^ b e d n N n ( a * bAab ~ 2^ a * bab,A ^ + Y, bed n N h ^ b/LN ^ abAa b ~ l^ abCl ^ A] +y (2J4) 

The Quantum Zero Range Model. Since no exclusion principle is operative for this 
model, we take its particles to be bosons. Thus, the operators a x and a* are unbounded 
here. In order to keep the formulation of the model in terms, exclusively, of bounded ones, 
we introduce the operators 

a x = (I + fU" 1 / 2 ^, a* = al(I + n,)" 1 / 2 (2.15) 

and note that, by Eqs. (2.6), (2.7) and (2.15), their actions on ifj(n) are given by the 
formula 

a x ^(n) = (1 - ^, )V(n x '-); a^(n) = iP(n x >+). (2.16) 
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Hence, a* and a x serve as bounded creation and annihilation operators. 

In order to formulate the quantum version of the generator of the dynamical semigroup 
of the model, we proceed along the same lines as for the simple exclusion model, simply 
replacing a x by a x . Thus we obtain the following formula for the quantum version of Eq. 
(2.3), as applied to bosons. 

QA = V' (g(h x ) 1/2 a x a y Aa y a x g(h x ) 1/2 - -\g(n x ) 1/2 aZaya y alg(nl /2 ),A] + ) + 

«o n(^(^6) 1/2 a^a 6 ^(n 6 ) 1 / 2 -^[^(n 6 ) 1 / 2 a^(n 6 ) 1/2 ,^] + ) + 

J2 bedi , N HWLN)(a b Aat-^[a b atA} + ). (2.17) 

2.3. The Quantum System £ as an Extension of the Classical one E c ; 

Since the subalgebra B of A is identified with that of the observables of E c ^, the 
following proposition establishes that, for the simple exclusion and zero range models, the 
dynamics of E induces an autonomous subdynamics on the classical observables B that is 
precisely that of the system E c j. In other words the quantum system E is an extension of 
the classical one, E c ;. 

Proposition 2.2. For the models under consideration, 

(1) The algebra B is stable under the semigroup 0(R + ); and 

(2) the restriction o/0(R + ) to B is just the dynamical semigroup </> c z(R + ) o/E c /. 

Proof . Since Q and Q c i are the generators of 4>(Ti + ) and c z(R + ), respectively, it 
suffices to show that , for the models concerned, Q c \ is just the restriction of Q to B. By 
Prop. 2.1, this condition is just that 

[gF(n)]iP(n) = [g c iF(n)]iP(n) V FeC(T), neT. (2.18) 

It is now a straightforward matter to check that, by Eqs. (2.5)-(2.8), (2.12), (2.14), (2.16) 
and (2.17), this condition is satisfied by both the simple exclusion and zero range models. 

Proposition 2.3. The dynamical transformations 0(R+) of the models under consid- 
eration commute with the gauge automorphisms 7(6*). Hence the dynamics of these models 
are locally gauge covariant. 

Proof. Since Q is the generator of 0(R + ), it suffices to show that Q commutes with 
7(0); and it follows from Eqs. (2.4)-(2.6), (2.10), (2.11), (2.15) and (2.17) that it does so. 

2.4. Steady States of E c; and E 

Assume now that E cZ has a unique steady state u) c \, as has been established for the 
simple exclusion and zero range models [1-5]. We shall now show that u c i extends to 
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a locally gauge invariant stationary state w of S. To this end we introduce the condi- 
tional expectation, P, of A onto B that defines PA as the mean over all the local gauge 
transformations 7(6*), i.e. 



r.27T 

PA = ' 



n xe n iV (27r)- 1 / d6 x 7 (0)A V A eA. (2.19) 
L Jo J 

We then define uj to be the state of E given by the formula 

uj(A) = u c i(PA) V AeA. (2.20) 

In view of Props. 2.2 and 2.3, it follows immediately from the last two equations that uj 
is indeed a locally gauge invariant stationary state of E, and that it is the only one that 
reduces to uj c i on £>. Moreover, in the case of the models under consideration, it is the 
only stationary state* of E, for the following reasons. Frigerio [12, Theorem 3.2] has shown 
that a quantum dynamical semigroup 0(R+) cannot admit more than one stationary state 
if the commutant of the operators H and {Vj} appearing in the formula (2.13) for its 
generator consists of the scalar multiples of the identity; and it follows easily from Eqs. 
(2.14) and (2.17) that, in view of the assumed strict positivity of the functions h and g, 
this condition is satisfied by the models under consideration. We remark that the steady 
states uj and uj c \ are ones of equilibrium or nonequilibrium according to whether or not the 
function h is constant over dQ^. 

In all cases, it follows from the above considerations that the quantum dynamical 
system E, as represented now by (^4, 0, uj), is an extension of the classical one (£>, <fi c i, uj c i). 

Since we shall be concerned with properties of the model in certain limits where N 
tends to infinity, we shall henceforth indicate the iV-dependence of E, E c ;, uj, uj c i, </>, <fi c i, Q 
and Qci by attaching the superscript (N) to these symbols. 

3. The Hydrodynamic Picture. 

We shall now investigate the large scale dynamical properties of the field n, with the 
aim of showing that it exhibits good hydrodynamical behaviour. To this end, we make the 
following two observations. 

(a) Since the field n of E^ is built from the observables observables n x , which are affiliated 
to the algebra £>, it follows from Props. 2.1 and 2.2 that the dynamics of this field reduces to 
that of its classical counterpart, n (of Section 2.1), as governed by the dynamical semigroup 
0if>(R + )ofE<r> 

(b) As remarked in the Note following Eq. (2.3), the evolution of the field n is diffusive. 
Hence, for a hydrodynamical description of this field on a length scale whose unit is L^r, 
the natural unit of the corresponding time scale is L 2 N . 

In view of these observations, we formulate the hydrodynamical picture of this field on 
macroscopic length and time scales whose units are Ln and L 2 N , respectively, as in Refs. 



In fact, the theory that follows does not depend on this uniqueness. 
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[1-5, 7, 8]. Thus, in this scaling, the dynamics of n is represented formally by the classical 
field 

qi N \x) = E o <t> {N HL 2 N t)n y 5(y - L N x) V xgO. yeK+. 

is therefore a £>'(0)-class distribution, in the sense of L. Schwartz [13]. To be precise, 
it is a continuous linear functional on the Schwartz space V(Q) of infinitely differentiate 
functions on R d with support in O; and its action on the latter space is given by the 
formula 



q ( t N \f) = L~ N d J2 pn ^(L^nyfiL-'y) V feV(Cl). (3.1) 

IN) 



In fact this field q\ does indeed exhibit good hydrodynamical properties since (cf. 
[2]), for appropriate initial states of E^i , the expectation value of q[ (/) converges 
in probability, as iV-^oo, to the smeared form f Q dxq t (x)f(x) of a smooth field which 
evolves according to a phenomenological equation of the form 

^ = A*(ft), (3.2) 

where the function $ is smooth and non-negative: in the case of the simple exclusion model 
it is the identity function. The spatial boundary condition for this evolution is given by 
the formula 

<$>(q t (x)) = h{x) V xedVL, teR+, (3.3) 

where h is the function that governs the boundary term in Eqs. (2.2) and (2.3). We denote 
by q the stationary solution of Eqs. (3.2) and (3.3). Evidently it is just the expectation 
value of ql for the nonequilibrium steady state u c i in the limit N^oo. 

We note that it follows from Eq. (3.2) that the linearised equation of motion for a 
small perturbations 5qt of q takes the form 

j t 5q t = C5q t (3.4) 

where 

£ = A[<&' (?(*))(.))] • (3.5) 

Hence, assuming that C is the generator of a one-parameter semigroup {Ti|tGR+} of linear 
transformations of X>'(0), the solution of Eqs. (3.4) is simply 

6q t = T t _ t0 6q t0 Vt>t (3.6) 



4. The Hydrodynamic Fluctuation Process 

We define the field which represent the fluctuations of q[ about its mean for 

the steady state , by the formula 

i{ N) if) = iV 1/2 bi N) (f) - U N) (/))] V fev(fl) . (4.1) 
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£( N ^ is thus a classical stochastic process for the state uj c i, indexed by Rx£>(fi). Our aim 
now is to verify that it satisfies the following conditions, which were the hypotheses on 
which the macrostatistical theory of [7,8] was based. 

(0) The hydrodynamic limit hypothesis. This asserts that the process converges 
in law to a stationary stochastic process £ as N^oo. 

(1) . The Regression Hypothesis. This asserts that the dynamical law governing devi- 
ations of the hydro dynamical variable qt from its steady state value is the same whether 
they arise from spontaneous fluctuations or from weak external perturbations. Thus, in 
view of the formula (3.6) for the perturbed hydrodynamics, the regression hypothesis is 
that 

£(6|6 ) = T t _ t0 £ t0 V t>to. (4.2) 
where E(.\£t ) denotes the conditional expectation, given 6 - 

(2) The Chaoticity Hypothesis. In order to specify this hypothesis, we introduce 
Nelson's [14] forward time derivative of 6, namely 

£>& = lim^or-^^+r - 616); (4-3) 

and we infer from Eqs. (4.2) and (4.3) that, since C is the generator of T(R + ), 

Dit = C£t- (4-4) 

By the definition (4.3), D^ t is the instantaneous expectation value of the rate of change of 
6- Accordingly, we designate f duD^ u to be the secular part of the increment (6 — 6) i n 
6 over the time interval [s, t]. Correspondingly, we designate the stochastic part of (6 — 6) 
to be the remaining part, wt yS , of this increment; and, in view of Eq. (4.4), this takes the 
form 

w t ,s = 6 - 6 - / duC^u- (4.5) 

J s 

Thus w is a process indexed by R^_xO. Our chaoticity hypothesis, which is designed to 
represent the stochasticity of this process, is that it is Gaussian and that its space-time 
correlations are of zero range, corresponding to ones of finite range on the microscopic 
scale. Thus the hypothesis is that w is Gaussian and that 

E(w t ,s(f) w t'A9)) = if either [s,t]f][s' ,t'} = or supp(/)nsupp((7) = 0. (4.6) 

( 3) The Local Equilibrium Hypothesis. We formulate the local properties of the process 
£ in terms of the transformation /— >/x ,e of V{Vt) defined by the formula 

f x ^ = e- d ' 2 f(e-\x-x )). (4.7) 

This transformation corresponds to the spatial rescaling by the factor e around the point 
xq. Further, in thermal equilibrium, the static two-point function for £ enjoy the properties 
[2,3] 

E{am(9)) = x(q) [ dxf(x)g(x), 
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where x represents the compressibility of the system; and 

E(ttC*(f)t( S ))=E(tif)tt£*g))= X m'(Q) [ dxV f{x).V g{x), 

where £* is the dual of C. By Eq. (4.6), these last two equations are equivalent to the 
following ones. 

E(H(fx M9xo,e))=x(<i) [ dxf(x)g(x), 

and 

e 2 E(Z(C*f X0 , £ )Z(g X0 , £ )) =e 2 E(((/ IO , e )((£\, £ )) =x(qW(q) [ dxV f(x).Vg(x). 

Jn 

The local equilibrium conditions, for fluctuations about nonequilibrium steady states, are 
just the limiting forms of these equations, as e decreases to zero, with q replaced by q(xo). 
Thus they are given by the formulae 

\im el0 E(Z(f X0 , e )ag X0 ,e))=x{<iM) [ dxf(x)g(x) (4.8) 

Jn 

and 

]h* eio( ?E(Z(f XOte )Z(£*g Xo , e )) =hm £i0 6 2 £;(e(£Vxo,e)e(^o ) e)) = 

X(q{x ))&(q{x )) [ dxVf{x).Vg{x). (4.9) 
Jn 

Note. These conditions represent local equilibrium on the hydrodynamic scale and are thus 
different from those formulated on the microscopic scale in Refs. [1, 2]. 

5. Verification of the Hypotheses 0-3 

We shall now show that the above hypotheses are verified by the simple exclusion and 
zero range models. 

The Simple Exclusion Model. The fluctuation process for this model was worked out 
in detail by Spohn [3] for the case where Q is the slab*, (0, l)xR d_1 . Here we shall confine 
our attention to the one- dimensional case, where O is the linear segment (0, 1). 

For this case the following results have been established [3] . 

(i) The function $ appearing in the phenomenological law (3.2) is just the identity function 
and correspondingly, by Eq. (3.5), the generator C is the Laplacian A, with Dirichlet 
boundary conditions. 

* Of course, for d > 1, the slab does not meet our condition that O be bounded. How- 
ever, it is a straightforward matter to extend our treatment and results to that situation. 



11 



(ii) The process ^ ' converges in law, as iV^oo, to a limit £ t that is governed by a 
Langevin equation 

dx t = C^ t dt + dw t , (5-1) 
where Wt is the Wiener process for which 

E([w t (f)-w s (f)][w t ,(g)-w s ,(g)}) = 

2 / dx X (q(x))^(q(x))Vf(x).Vg(x)\[s,t} n [s',t'}\ V f,geV(Q), t,s(<t),t,s'(<t)eIL + 

(5.2) 

and 

E([w t - w s ]\Q = for t>s>u.. (5.3) 

(iii) The static two-point function takes the form 

E(am(9)) = f dx X (q(x))f(x)g(x) + [h(l) - h(0)} f dxf(x)A- 1 g(x) 1 (5.4) 

JO JO 

where 

x(q) = q( 1 -q) (5-5) 

and h is the function appearing in the boundary term of Eq. (2.2). 

The result (ii) immediately substantiates the hypotheses (0) and (2). Moreover, since 
C is the generator of £(R + ), it follows from Eq. (5.1) that 

6 = T t _ t0 & + / T t - U dw{u) V t>t (5.5) 

and hence, by Eq. (5.3), that the hypothesis (1) is also fulfilled. Finally, the local equi- 
librium properties (4.8) and (4.9) are simple consequences of the formulae (4.7) and (5.4), 
which signifies that the model also satifies hypothesis (3). 

The Zero Range Model. A key property of this model is that its steady state takes 
the simple product form [15, 5] 

^IP = ®x£n N m Xil}{x ), (5.6) 

where m x ^ x ^ is a probability measure on the functions of n x that depends on the value 
of the stationary field q at the site x. 

It follows now from a straightforward adaptation of the argument* of [2, Ch. 11] that, 
for this model too, the process converges in law to a Gaussian process £, which is 

* For that argument, as applied to the present situation, £ f lies in the Sobolev space 
H- r {Q) := {/ : O^R| / n dxf(l - A)~ n f < oo} for sufficiently large r(eN). This space is 

1/2 

a subspace of V(Q). The operators £* and [$(g(x))] V on H- r play the roles of those 
denoted, in Gothic script, by A and £>, respectively, in [2]. 
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also represented by Eqs. (5.1)-(5.3) , though with C now given by Eq. (3.5). Hence, by 
the same argument as for the simple exclusion model, we see that this model satisfies the 
hypotheses (0), (1) and (2). Further, it follows from (0) that 

and hence, by Eqs. (4.7) and (5.6), that the conditions (4.8) and (4.9) are fulfilled, with 
x{q( x )) the variance of the particle number at the site Lnx in the single particle state 

6. Concluding Remarks 

In this article we have constructed a quantum stochastic model that fulfills the hy- 
potheses of our general macrostatistical picture [7,8] of nonequilibrium steady states. A 
rather unphysical feature of this model is that the density field q[ N ^ is classical not only at 
the hydrodynamical level but also at the microscopic one. Thus the problem of construct- 
ing a quantum model whose classical properties emerge only at the hydrodynamical and 
thermodynamical levels remains a challenging and interesting one. 
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